Discrimination of supraventricular tachycardia and ventricular tachycardia events

ABSTRACT

A method and system for discrimination of supraventricular tachycardia and ventricular tachycardia events. Morphological features points are extracted from normal sinus rhythm (NSR) complexes and used to generate a NSR template. A numerical convolution is performed using the NSR template and the feature points for each sensed NSR to give a NSR filter output. Using a plurality of NSR complexes, a median NSR filter output template is determined, where the median NSR filter output template has a median value for each value in the NSR filter output. The median NSR filter output template is then used during a tachycardia event to distinguish tachycardia events as either ventricular tachycardia events or supraventricular tachycardia events.

FIELD OF THE INVENTION

This invention relates generally to medical devices, and more particulary to a system and method for discriminating supraventricular tachycardia from ventricular tachycardia during a tachycardia event.

BACKGROUND

The heart is generally divided into four chambers, the left and right a trial chambers and the left and right ventricular chambers. As the heart beats, the atrial chambers and the ventricular chambers go through a cardiac cycle. The cardiac cycle consists of one complete sequence of contraction and relaxation of the chambers of the heart. The terms systole and diastole are used to describe the contraction and relaxation phases the chambers of the heart experience during a cardiac cycle. In systole, the ventricular muscle cells contract to pump blood through the circulatory system. During diastole, the ventricular muscle cells relax, causing blood from the atrial chamber to fill the ventricular chamber. After the period of diastolic filling, the systolic phase of a new cardiac cycle is initiated.

Through the cardiac cycle, the heart pumps blood through the circulatory system. Effective pumping of the heart depends upon five basic requirements. First, the contractions of cardiac muscle must occur at regular intervals and be synchronized. Second, the valves separating the chambers of the heart must fully open as blood passes through the chambers. Third, the valves must not leak. Fourth, the contraction of the cardiac muscle must be forceful. Fifth, the ventricles must fill adequately during diastole.

When the contractions of the heart are not occurring at regular intervals or are unsynchronized the heart is said to be arrhythmic. During an arrhythmia, the heart's ability to effectively and efficiently pump blood is compromised. Many different types of arrhythmias have been identified. Arrhythmias can occur in either the atrial chambers or in the ventricular chambers of the heart.

Ventricular fibrillation is an arrhythmia that occurs in the ventricle chambers of the heart. In ventricular fibrillation, various areas of the ventricle are excited and contract asynchronously. During ventricular fibrillation the heart fails to pump blood. Since no blood is pumped during ventricular fibrillation, the situation is fatal unless quickly corrected by cardiac conversion. Ventricular tachycardia is another arrhythmia that occurs in the ventricular chambers of the heart. Ventricular tachycardia is a very serious condition. Ventricular tachycardias are typified by ventricular rates between 120-250 and are caused by disturbances (electrical or mechanical) within the ventricles of the heart. During a ventricular tachycardia, the diastolic filling time is reduced and the ventricular contractions are less synchronized and therefore less effective than normal. Ventricular tachycardias must be treated quickly in order to prevent the tachycardia from degrading into a life threatening ventricular fibrillation.

Arrhythmias that occur in the atrial chambers of the heart are referred to generally as supraventricular tachycardias. Supraventricular tachycardias include atrial tachycardias, atrial flutter and atrial fibrillation. During certain supraventricular tachycardias, aberrant cardiac signals from the atria drive the ventricles at a very rapid rate. Such a situation occurs during paroxysmal atrial tachycardia. This condition begins abruptly, lasts for a few minutes to a few hours, and then, just as abruptly, disappears and the heart rate reverts back to normal.

Cardioverter-defibrillators, such as implantable cardioverter-defibrillators (ICDs), have been shown to be effective in reducing the incidence of sudden cardiac death. Sudden cardiac death is typically caused by either ventricular tachycardia or ventricular fibrillation. Cardioverter-defibrillator systems operate by sensing and analyzing cardiac signals and applying electrical energy to the heart when either a ventricular tachycardia or ventricular fibrillation is detected. One common way cardioverter-defibrillators detect cardiac arrhythmias is to sense and analyze the rate of ventricular contractions. When the ventricular rate exceeds a programmed threshold value, the cardioverter-defibrillator applies electrical energy in one or more specific patterns to treat either the ventricular tachycardia or ventricular fibrillation.

Rapid ventricular rhythms, however, can occur in the presence of a supraventricular tachycardia. As previously mentioned, one example is during paroxysmal atrial tachycardia. In this situation, treating the ventricles with electrical energy is inappropriate as the treatment does not address the precipitating factor of the rapid ventricular rate. Therefore, a need exists for reliably assessing and determining the origin of a rapid ventricular rate. By reliably discriminating the origin of the rapid ventricular rate, more appropriate and effective therapies can be applied to treat the heart.

SUMMARY OF THE INVENTION

The present subject matter discloses a method and a system for discriminating, or classifying supraventricular tachycardias (SVT) from malignant ventricular tachycardias (VT) during a tachycardia event. In one embodiment, the present subject matter is implemented in an implantable cardioverter defibrillator. By using the method of the present subject matter, the implantable defibrillator assesses and determines the origin of a rapid ventricular rate, allowing the implantable defibrillator to reduce the number of inappropriate therapies delivered to the heart.

In one embodiment, QRS-complexes are sensed during normal sinus rhythm (NSR). A plurality of feature points are located on the sensed NSR QRS-complexes based on morphological features of the individual NSR QRS-complexes. The plurality of feature points from the NSR QRS-complexes are then used to determine a NSR template. In one embodiment, a plurality of NSR QRS-complexes are used to determine the NSR template.

In one embodiment, the NSR template includes a median value for each of the plurality of feature points taken along the NSR QRS-complex. A numerical convolution is then preformed on the values of the NSR template. A numerical convolution is also preformed on the plurality of feature points for each of the plurality of the NSR QRS-complexes. This process gives a NSR filter output for each of the NSR QRS-complexes. Using the NSR filter output for each NSR QRS-complex, a median NSR filter output template is determined. In one embodiment, the median NSR filter output template includes the median values of the NSR filter output values for each NSR QRS-complex.

Once the median NSR filter output template has been determined, the system senses for the occurrence of a tachycardia event. When a tachycardia event is detected, the system senses the tachycardia complexes. In one embodiment, QRS-complexes are extracted, or sampled, from the tachycardia complexes in the sensed signals. The plurality of feature points are then located in the QRS-complexes. The feature points located in the QRS-complexes during the tachycardia event are based on morphological features of the QRS-complex. In one embodiment, the morphological features taken from the QRS-complex during the tachycardia episode are from the same relative position as the morphological features taken along the NSR QRS-complex.

A tachycardia complex output is then determined by performing a numerical convolution of the median NSR filter output template with the plurality of feature points from a QRS-complex of a tachycardia complex sensed during the tachycardia event. The differences between the values of the tachycardia complex output and the median NSR filter output template are summed to give a sum of residual value. The sum of residual (SOR) value is then compared to a sum of residual (SOR) threshold value, and when the SOR value is greater than or equal to the SOR threshold value the tachycardia complex is classified as a ventricular tachycardia complex. When the SOR value is less than the SOR threshold value the tachycardia complex is classified as a supraventricular tachycardia complex. When the number of tachycardia complexes classified as either ventricular tachycardia complexes or supraventricular tachycardia complexes exceeds a predetermined threshold value, the tachycardia event is classified as the tachycardia event that exceeded the threshold value (i.e., as either a ventricular tachycardia or a supraventricular tachycardia).

These and other features and advantages of the invention will become apparent from the following description of the preferred embodiments of the invention.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic of an implantable medical device;

FIG. 2 is a flow chart illustrating one embodiment of the present subject matter;

FIG. 3 is a flow chart illustrating one embodiment of the present subject matter;

FIG. 4 shows one embodiment of an electrocardiogram of a normal sinus rhythm complex;

FIG. 5 is a flow chart illustrating one embodiment of the present subject matter;

FIGS. 6A, 6B and 6C show plots of sensed cardiac complexes as a function of time for three cardiac conditions;

FIGS. 6D, 6E and 6F show plots of SOR values as a function of sensed cardiac complex for the sensed cardiac complexes of FIGS. 6A, 6B and 6C; and

FIG. 7 is one embodiment of a block diagram of an implantable medical device according to the present subject matter.

DETAILED DESCRIPTION

In the following detailed description, reference is made to the accompanying drawings which form a part hereof and in which is shown by way of illustration specific embodiments in which the invention can be practiced. These embodiments are described in sufficient detail to enable those skilled in the art to practice and use the invention, and it is to be understood that other embodiments may be utilized and that electrical, logical, and structural changes may be made without departing from the spirit and scope of the present invention. The following detailed description is, therefore, not to be taken in a limiting sense and the scope of the present invention is defined by the appended claims and their equivalents.

The embodiments illustrated herein are demonstrated in an implantable cardiac defibrillator (ICD), which may include numerous defibrillation, pacing, and pulse generating modes known in the art. However, these embodiments are illustrative of some of the applications of the present system, and are not intended in an exhaustive or exclusive sense. The concepts described herein can be used in a variety of applications which will be readily appreciated by those skilled in the art upon reading and understanding this description. For example, the present system is suitable for implementation in a variety of implantable and external devices.

As discussed above, ICDs can deliver inappropriate therapy to the ventricles of the heart. One example is when a supraventricular tachycardia causes a rapid ventricular rate. The ICD senses the rapid ventricular rate and responds by treating the heart (i.e., delivering electrical energy to the ventricles) as if a ventricular tachycardia were occurring. However, when the cause of the rapid ventricular rate is a supraventricular tachycardia treating the ventricles is inappropriate as the therapy is not directed at the source of the arrhythmia. Looking at it another way, the ICD is attempting to treat the symptoms, not the cause of the rapid ventricular rate.

The present subject matter allows for the origin, or the source, of a rapid ventricular rate to be distinguished. In one embodiment, the present subject matter distinguishes the origin of a ventricular tachycardia as either the result of a supraventricular tachycardia (i.e., an arrhythmia in the atria is driving the ventricles at a rapid rate) or the result of a disturbance in the ventricles (i.e., a ventricular tachycardia).

Rapid ventricular rates can include ventricular rates that are between 100-250 beats per minute, as are typical during ventricular tachycardia. Ventricular rates of greater than 250 beats per minute are also considered to be within the present subject matter.

The present subject matter utilizes sensed cardiac complexes in determining, or distinguishing, the source of rapid ventricular rates. Cardiac complexes include the electrical excitatory impulses, or action potentials, that are sensed from the heart is it goes through the cardiac cycle. In one embodiment, the sensed cardiac complexes are electrocardiogram signals sensed from the beating heart. Major features of an electrocardiogram signal include the P-wave, the QRS-complex and the T-wave which are caused by the atrial depolarization, the ventricular depolarization and the ventricular repolarization, respectively.

In one embodiment, cardiac complexes are sensed and analyzed using an ICD system. In one embodiment, the ICD system uses a transvenous lead system to allow sensing of the cardiac action potentials. The transvenous lead system can include a rate-sensing electrode and at least one defibrillation electrode positioned on the transvenous lead. Cardiac action potentials sensed using defibrillation electrodes are typically referred to as far-field signals (or morphology signals) and cardiac action potentials sensed using rate sensing, or pacing, electrodes are typically referred to as near-field signals (or rate signals).

In one embodiment, the rate-sensing electrode is a pacing tip electrode positioned at the distal end of the transvenous lead system. Other types of rate sensing electrodes are also considered appropriate to use with the present subject matter. Examples of other types of rate sensing electrodes include ring electrodes, both annular and semi-annular, as are known in the art. Rate sensing using the transvenous lead system can also be accomplished either through unipolar or bipolar sensing methods, as are known.

In one embodiment, the transvenous lead system can have a single defibrillation electrode. When a single defibrillation electrode is present, the ICD uses unipolar sensing and defibrillation, as is known. In one embodiment, a transvenous catheter with a single defibrillation electrode is implanted with the rate sensing electrode and the defibrillation electrode positioned within the right ventricular chamber of the heart. In an alternative embodiment, the transvenous lead can have two, or more, defibrillation electrodes. When two defibrillation electrodes are present, the ICD system can preform bipolar sensing the cardiac action potentials. In bipolar sensing, cardiac action potentials are sensed between the two defibrillation electrodes, and defibrillation electrical shocks are delivered between the two defibrillation electrodes. The ICD system of the present subject matter can also use endocardial and patch electrodes as are known.

In one embodiment, the ICD employs an single body lead catheter sold under the trademark ENDOTAK (Cardiac Pacemaker, Inc./ Guidant Corporation, St. Paul, Minn. having a pacing tip electrode and two defibrillation coil electrodes. One example of such a system is shown in FIG. 1. ICD 100 is coupled to catheter 110, which is implanted to receive signals from heart 120. The catheter 110 also may be used for transmission of pacing and/or defibrillation signals to the heart 120. In an alternative embodiment, a three defibrillation electrode system is employed, wherein the housing of the implantable system is used as a third defibrillation electrode.

In one example, the ICD 100 senses cardiac signals from the heart 120. When the ICD 100 detects the occurrence of an arrhythmic event, the ICD 100 analyzes the sensed arrhythmic complexes (i.e., the cardiac signals) of thc arrhythmic event. During the analysis, the ICD compares the sensed arrhythmic complexes to cardiac signals sensed and recorded during the patient's normal sinus rhythm. Based on the comparison, the ICD 100 is able to distinguish SVT events from VT events and, depending upon the ICD's classification of the arrhythmic event, to provide appropriate therapy to treat the heart 120.

FIG. 2 shows another embodiment of a system for distinguishing SVT events from VT events during a tachycardia event. At 200, the system senses cardiac signals representative of electrical cardiac activity. At 210, the system analyzes the sensed cardiac signals to detect the occurrence of a tachycardia event in the heart. In one embodiment, the tachycardia event are ventricular arrhythmic events.

When tachycardia event is detected at 210, the system continues to sense the tachycardia complexes. At 220, the sensed tachycardia complexes are then compared to a template in order to determine the classification of the sensed tachycardia complex. In one embodiment, the template is determined from cardiac complexes sensed during normal sinus rhythm. For example, features from the cardiac complexes sensed during normal sinus rhythm are used to create a normal sinus rhythm template. Features from the tachycardia complexes sensed during the tachycardia event are then compared to the normal sinus rhythm template. In one embodiment, the features from the sensed tachycardia complexes correspond to the features taken from the cardiac complexes sensed during normal sinus rhythm. In other words, features taken from the tachycardia complexes are from the same relative location as those take from the normal sinus rhythm cardiac complexes.

In one embodiment, the comparison of the sensed tachycardiac complexes and the normal sinus rhythm template at 220 determines how similar or dissimilar the sensed tachycardiac complexes are from the patient's normal sinus rhythm complexes. Based on the comparison, the system at 230 then classifies the sensed tachycardiac complexes as either ventricular tachycardiac signals or supraventricular tachycardia signals. As the tachycardiac complexes are classified, the system counts the occurrence and the classification of each tachycardiac complex. At 240, the system then determines whether the number of sensed ventricular tachycardia complexes has exceeded a predetermined threshold value. If the threshold has been exceeded, the system proceeds to 250 where a ventricular tachycardia event is declared and therapy to treat the ventricular tachycardia is delivered to the patient's heart. If the threshold has not been exceeded, the system returns to 210, where additional cardiac complexes are sensed and classified by the system.

FIG. 3 shows an embodiment for determining a normal sinus rhythm (NSR) template. At 300, the system senses cardiac complexes during normal sinus rhythm (NSR). In one embodiment, the NSR complexes are sensed and recorded using the ICD under the supervision of the patient's attending physician. In one embodiment, the system samples, or senses, QRS-complexes of the NSR complexes.

As each QRS-complex is sensed it is isolated, or windowed, for analysis. At 310, a plurality of feature points are located on the QRS-complex. In one embodiment, the plurality of feature points are morphological features of the QRS-complex. For example, morphological features extracted from the QRS-complexes include the amplitude values of peaks and valleys (or maxima and minima) acquired by a process called feature extraction. During feature extraction, each NSR complex is isolated according to a known morphological template. In one embodiment, the morphological template operates to detect the activation of an heart beat (such as the occurrence of an R-wave), at which point the electronic control circuitry of the implantable medical device analyzes the complex associated with the signal indicating the activation of the heart beat. In one embodiment, a threshold value or a detection criterion is used to indicate the activation of the heart beat. Once a heart beat has been detected, the feature extraction derives the plurality of feature points from the morphological features of the NSR complex. The extracted feature points for each NSR complex create a vector. In one embodiment, the vector includes a set of numbers, where each number is associated with a particular morphological point along the sensed NSR complex. Other types of features known in the art can also be extracted and used in developing the NSR template.

In one embodiment, two or more morphological features are extracted and used from the NSR complexes in developing the NSR template. For example, four features can be extracted from each NSR complex in developing the NSR template. In one embodiment, the features having a numerical value which are characteristic of the morphological position along the cardiac complex. So, at 310, as the NSR complexes are sensed, a set of four features are extracted, and a four feature NSR complex is created and stored for each of the NSR complexes. FIG. 4 shows one embodiment of a NSR complex 400. The NSR complex 400 is processed to determine the amplitudes of peaks 402 and valleys 404 in the QRS-complex 406. In one embodiment, the peaks 402 and valleys 404 are determined by determining major inflection points in the QRS-complex as represented in FIG. 4. The resulting values of the peaks 402 and valleys 404 provides a four dimensional NSR complex vector, [1A, 1B, 1C, 1D], where the number “1” represents the number of the sensed NSR complex. As previously mentioned, other features from cardiac complexes known in the art can be used in developing vectors. By way of example, other features can include the start or ending of a cardiac complex as detected by a predetermined deviation from a baseline signal or by detecting a predetermined decrease or increase in the slope of the cardiac signal. Thus, the present subject matter is in no way limited to vectors developed from maximum and minimum deflection points along signals of cardiac complexes.

Referring again to FIG. 3, at 320, a NSR template is determined from a plurality of the NSR complex vectors. In one embodiment, the NSR template is a median value computed from the corresponding NSR complex vectors. For example, the NSR template [TA, TB, TC, TD] is determined from a series of NSR complex vectors [1A, 1B, 1C, 1D], [2A, 2B, 2C, 2D], [3A, 3B, 3C, 3D]. . . , [nA, nB, nC, nD], where TA equals the median value of (1A, 2A, 3A, . . . nA), TB equals the median value of (1B, 2B, 3B, . . . nB), TC equals the median value of (1C, 2C, 3C, . . . nC), and TD equals the median value of (1D, 2D, 3D, . . . nD). In an alternative embodiment, the NSR template feature value is an average value computed from the corresponding NSR complex vectors. The number of NSR complexes used in determining TA, TB, TC, TD is a programmable value. In one embodiment, the number of NSR complexes used in determining the NSR template is a programmable value in the range of 10 to 100 NSR complexes.

At 330, the NSR template is used to create the matched filter impulse response (h(t)). In one embodiment, the matched filter impulse response is defined to be:

h(t)=[TD, TC, TB, TA]

The h(t) is “matched” to the NSR template. A matched filter is a specific type of filter designed to maximize the output signal-to-noise ratio. The matched filter effectively correlates an input signal with a stored replica of a signal of interest. The impulse response of the matched filter is the signal of interest time reversed and possibly time shifted. When a signal of interest is detected by the matched filter, the output signal-to-noise ratio will be maximized.

In one embodiment, the h(t) is used to filter features from tachycardia complexes sensed during a tachycardia event. In an additional embodiment, the system at 340 performs a numerical convolution of the NSR template and the NSR complexes to give a NSR filter output. The system at 350 then determines a NSR filter output template from the NSR filter output. In one embodiment, features are extracted from normal sinus rhythm complexes to create the NSR filter output template. In one embodiment, the features extracted from the NSR complexes create a four element vector as previously discussed. The extracted features are taken from the same relative position along the tachycardia complexes as the features taken from the NSR complexes. After extracting the features from the sensed tachycardia complexes, the feature values are stored for processing or filtering.

In one embodiment, as tachycardia complexes are sensed during a tachycardia event, the system extracts the values of the four features along the tachycardia complexes. In one embodiment, the values of the four features are referred to as vectors. As the vector is determined for each sensed tachycardia complex, the system filters the complex with respect to the four features of the normal sinus rhythm template. The result is a seven (7) element filter output for each of the sensed cardiac complexes.

In one embodiment, the vector for each of the NSR cardiac complexes are convolved with h(t). In one embodiment, this procedure produces the seven (7) element filter output for each of the sensed cardiac complexes. The seven element filter output values are then used to determine the NSR filter output template. In one embodiment, the NSR filter output template is a median value computed from a plurality of normal sinus rhythm complexes, where the median NSR filter output template has a median value for each value in the NSR filter output. For example, the median NSR filter output template [TE, TF, TG, TH, TI, TJ, TK] is determined from aplurality of NSR complex vectors [1E, 1F, 1G, 1H, 1I, 1J, 1K], [2E, 2F, 2G, 2H, 2I, 2J, 2K], [3E, 3F, 3G, 3H, 3I, 3J, 3K]. . . , [nF, nG, nH, nI, nJ, nK], where TE equals the median value of (1E, 2E, 3E, . . . nE), TF equals the median value of (1F, 2F, 3F, . . . nF), TG equals the median value of (1G, 2G, 3G, . . . nG), TH equals the median value of(1H, 2H, 3H, . . . nH), TI equals the median value of (1I, 2I, 3I, . . . nI), TJ equals the median value of (1J, 2J, 3J, . . . nJ), and TK equals the median value of (1K, 2K, 3K, . . . nK). In an alternative embodiment, the NSR filter output template is an average value computed from the corresponding NSR complex vectors.

In one embodiment, the method for determining a normal sinus rhythm (NSR) template is preformed with an implantable cardioverter defibrillator. The ICD system discussed herein is one example of an appropriate system for determining a patient's NSR template. In an alternative embodiment, other implantable medical devices, such as implantable defibrillators, or external defibrillators can be used to determine the NSR template and implement the present subject matter.

FIG. 5 shows an embodiment of a system for distinguishing SVT events from VT events during a tachycardia event. At 500, cardiac signals representative of electrical cardiac activity are sensed. In one embodiment, the cardiac signals are sensed by an endocardial lead system of an ICD as previously discussed. The cardiac signals include cardiac complexes which are portions of the complete cardiac cycles. In one embodiment, the sensed cardiac complexes include the QRS-complex of the cardiac cycle. The system analyzes the sensed cardiac complexes to determine if a tachycardia event is occurring. In one embodiment, the system determines the occurrence of a tachycardiac event by analyzing the sensed cardiac rate. A cardiac rate that exceeds a predetermined threshold indicates the occurrence of a ventricular tachycardia. In one embodiment, the predetermined threshold is for cardiac rates of between 150-250 beats per minute. In an alternative embodiment, the predetermined threshold is a lower rate zone in an ICD that is configured with multiple rate-zones. Other methods of determining the occurrence of tachycardia event which are known in the art may be used without departing from the present system.

When a tachycardia event is detected, the system proceeds to 504. At 504, the system samples the tachycardia complexes during the tachycardia event. Feature points along the tachycardia complexes are located on the QRS-complex of the tachycardia complexes at 508. In one embodiment, the feature points located along the sensed tachycardia complexes are at the same relative position as feature points located along QRS-complexes sensed during normal sinus rhythm. For example, the QRS-complex is isolated, or windowed, for analysis and vector values derived for the complex as previously discussed.

In one embodiment, the set of features extracted from the QRS-complex of the cardiac complex are processed using the NSR template. In one embodiment, the set of features extracted from the QRS-complex are convolved with the matched filter impulse response (h(t)) of the NSR template at 512 to give a tachycardia complex output. In one embodiment, the h(t) used at 512 is derived from features extracted from a patient's normal sinus rhythm complexes.

In processing the tachycardiac complexes with h(t), the system extracts features from the morphology of the tachycardiac complexes. In one embodiment, the number of features extracted from the morphology signal of the tachycardiac complexes is equal to the number of features extracted from the morphology signals of the sensed normal sinus rhythm complexes. So, for example, four features are extracted from the morphology of the sensed normal sinus rhythm complexes and so four features are extracted from each of the sensed tachycardiac complexes. The extracted features are taken from the same relative position along the sensed tachycardiac complexes as those taken from the complexes sensed during normal sinus rhythm. After extracting the features from the sensed tachycardiac complexes, the feature values are stored for processing, or filtering (including the use of a matched filtering), with the normal sinus rhythm template at 512.

In one embodiment, as the tachycardiac complexes are sensed during the tachycardia event, the system measures values for each of the features along the tachycardiac complex. In one embodiment, the values of the features are referred to as vectors. As the vector values are determined for each sensed tachycardiac complex, the system filters the signal with respect to the features of the normal sinus rhythm template (or the matched filter impulse response) to give the tachycardia complex output. For example, when four features are sensed from the normal sinus rhythm complexes and the tachycardiac complexes, the system filters the tachycardiac complexes with respect to the four features of the normal sinus rhythm template (or matched filter impulse response). When four features are used, the result of the filtering process is a seven (7) element filter output for each of the sensed tachycardiac complexes. In this embodiment, the seven element filter output is the tachycardia complex output for the sensed cardiac complex.

In one embodiment, the filtering process is similar to a mathematical convolution of the normal sinus rhythm template with the features extracted from the sensed cardiac signals. In one embodiment, the process of numerical convolution can be envisioned as the interaction between two strips. In one embodiment, the numerical convolution of y(t)=f(t)*h(t) can be described where the sequence of f[1], f[2], f[3], . . . , [samples for f(t)] are written on the upper strip, and the sequence h[1], h[2], h[3], . . . , [samples for h(t)], are written on the lower strip. The upper strip is fixed, and the lower strip is folded (inverted) about t=1 and then moved from left to right, one slot at a time, with the f[1] and h[1] slots coinciding at t=1. To compute y(t), the lower strip is shifted by t slots and multiply the sample values of f(t) by the samples values of h(t) lying in the adjacent slots. The product of the multiplied samples are then summed together. This procedure is then repeated for each time the lower strip is shifted along the upper strip. For example, in filtering the four features of the cardiac complex for the present system, the result is a seven member vector (the tachycardia complex output for the sensed cardiac complex), where the seven member vector is determined from the calculation of:

[h(1)f(1), h(1)f(2)+h(2)f(1), h(1)f(3)+h(2)f(2)+h(3)f(1), h(1)f(4)+h(2)f(3)+h(3)f(2)+h(4)f(1), h(2)f(4)+h(3)f(3)+h(4)f(2), h(3)f(4)+h(4)f(3), h(4)f(4)]

At 516, the system then sums the difference between the tachycardia complex output and the median NSR filter output template. In one embodiment, the resulting difference is given the term sum of residuals (SOR) value. In one embodiment, a SOR value is calculated for each of the sensed tachycardiac complexes. One way to think of the SOR is as a differentiation value, where the SOR is used as a comparison of magnitudes at specific points along morphological signals. In one embodiment, the SOR is calculated from the seven-element output vectors from the sensed tachycardiac complexes and the normal sinus rhythm output template.

As the SOR can be thought of as the sum of differences between the tachycardia complex output and the median NSR filter output template, there are many ways in which the differences can be calculated. For example, the SOR is determined by taking the absolute value of the differences of corresponding elements from a seven-element output vector and the single, seven-element NSR output template. These quantities are then added together to produce a scalar quantity for each complex. This computation can also be described mathematically as: ${{\text{Sum of Residuals} = {\sum\limits_{i = 1}^{7}{{N_{i} - T_{i}}}}},\text{Where:}}\quad$ N = seven-elementNSR  output template. T = seven-element tachycardia  output vector.

In an additional embodiment, the SOR value can also be calculated by determining the sum of the square of the differences between N_(i) and T_(i). Alternatively, the SOR value can be calculated by summing the difference between N_(i) and T_(i) Other mathematical methods of quantifying differences between N_(i) and T_(i) exist and are considered within the scope of the present invention.

Referring now to FIGS. 6A, 6B and 6C there are shown sensed cardiac complexes plotted as a function of time. FIG. 6A shows cardiac complexes sensed during normal sinus rhythm; FIG. 6B shows cardiac complexes sensed during ventricular tachycardia; and FIG. 6B shows cardiac complexes sensed during atrial fibrillation.

FIGS. 6D, 6E and 6F show SOR values plotted as a function of sensed cardiac complexes. FIG. 6D shows SOR values calculated for the cardiac complexes plotted in FIG. 6A; FIG. 6E shows SOR values calculated for the cardiac complexes plotted in FIG. 6B; and FIG. 6F shows SOR values calculated for the cardiac complexes plotted in FIG. 6C.

Comparing the calculated SOR values for each of FIGS. 6D, 6E and 6F reveal differences in SOR values for the cardiac complexes sensed during the ventricular tachycardia event and the atrial fibrillation event as compared to cardiac complexes sensed during normal sinus rhythm. Additionally, differences in SOR values can be seen for the cardiac complexes sensed during the ventricular tachycardia event as compared to those of the atrial fibrillation event.

Referring again to FIG. 5, after the sum of the differences between the tachycardia complex output and the median NSR filter output template is computed at 516, the system compares the numerical difference to a threshold value at 524. In one embodiment, the SOR value for each tachycardiac complexes is compared to a SOR threshold value. In one embodiment, when the numerical difference is less than the threshold value, the cardiac complex is classified as a supraventricular tachycardia (SVT) complex at 528. In one embodiment, when the numerical difference is greater than, or equal to, the threshold value, the cardiac complex is classified as a ventricular tachycardia (VT) complex at 532.

As the sensed tachycardiac complexes are categorized, the system records the number of ventricular tachycardia complexes and SVT complexes that have been categorized during the tachycardia event at 536. In one embodiment, the predetermined threshold is an x out of the last y signals counter. In one embodiment, the values for x and y are programmable. In one embodiment, x has programmable integer values of greater than 3, where the values can be in the range of 3 to 10, where, in an additional embodiment, 5 is an acceptable value. In one embodiment, y has a programmable integer value of greater than 8, where the values can be in the range of 8 to 30, where 10 is an acceptable value. In an alternative embodiment, the system determines a percentage of VT complexes during the tachycardia event. When the percentage of the VT complexes exceeds a predetermined percentage threshold, the system declares the occurrence of a ventricular tachycardia. In one embodiment, the predetermined percentage threshold is a programmable value in the range of 30 to 100 percent, where 50 percent is an acceptable value.

When the number of VT complexes exceeds the predetermined threshold, a VT event is declared. A signal is then provided to the system (e.g., an implantable medical device) to deliver ventricular tachycardia therapy at 540 to a heart when a ventricular tachycardia episode is declared. In an alternative embodiment, if the number of SVT complexes exceeds the predetermined threshold, an SVT event is declared. In one embodiment, when an SVT event is declared, a signal is provided to the system (e.g., an implantable medical device) to deliver supraventricular tachycardia therapy to the heart. In an alternative embodiment, therapy is not delivered to the supraventricular region of the heart, but rather the system continues to monitor the cardiac condition and provides treatment only when a ventricular tachycardia is determined. If at 536 the number of classified cardiac complexes does not exceed the predetermined threshold, the system returns to 508 to sense and classify the next tachycardiac complex.

FIG. 7 shows one embodiment of an implantable cardiac defibrillator (ICD) 700, which may include numerous defibrillation, pacing, and pulse generating modes known in the art. An endocardial lead is physically and electrically coupled to the ICD 700. The endocardial lead can include at least one pacing electrode and at least one defibrillation coil electrode as are known. In one embodiment, the endocardial lead is an ENDOTAK lead as previously described.

FIG. 7 discloses ICD 700 which includes input circuitry 710. In one embodiment, input circuitry 710 includes a first amp 712 and a second amp 714. The first amp 712 receives rate-signals or near-field signals through the at least one pacing electrode. In one embodiment, the rate-signals are sensed using a unipolar configuration, where the cardiac signals are sensed between the at least one pacing electrode and the housing 716 of the ICD 700. Alternatively, bipolar sensing is accomplished between two or more pacing electrodes on one or more endocardial leads. The second amp 714 receives morphology-signals, or far-field signals, from at least two defibrillation coil electrodes located on the endocardial lead.

An R-wave detector 720 receives the rate-signals from the first amp 712. The R-wave detector 720 detects R-waves from the rate-signals being received by the first amp 712 and conveys information relating to the cardiac rate to a microprocessor 724 by a data bus 726. A morphology analyzer 730 receives morphology signals from the second amp 714. In one embodiment, the morphology analyzer 730 extracts a plurality of feature points from sensed cardiac complexes. A template generating circuit 734 is coupled to the signal morphology analyzing circuit 730 by the bus 726. The template generating circuit 734 receives the extracted plurality of feature points. In one embodiment, the template generating circuit 734 generates a normal sinus rhythm template from sensed normal sinus rhythm complexes.

A filter output response circuit 740 is coupled to the template generating circuit 734 by bus 726. The filter output response circuit 740 creates a normal sinus rhythm filter output by performing a numerical convolution on the NSR template and the plurality of feature points for each of the plurality of the NSR complexes. The filter output response circuit 740 also determines a median normal sinus rhythm output template from the plurality of normal sinus rhythm complexes, where the median NSR filter output template has a median value for each value in the NSR filter output.

In one embodiment, when, during a tachycardia event, the input circuitry 710 receive a QRS-wave signal from a tachycardia complex, the signal morphology analyzing circuit 730 locates the plurality of feature points on the QRS-complex based on morphological features of the QRS-complex. The filter output response circuit 740 then performs a numerical convolution of the NSR template with the plurality of feature points on the QRS-complex to give a tachycardia complex output. The filter output response circuit 740 then sums a numerical difference between the values of the tachycardia complex output and the median NSR filter output template.

The summed difference is then received by the microprocessor 724 where the calculated value is compared to the predetermined sum of residual threshold value. During a tachycardia episode, as the system senses each cardiac complex the ICD 700 classifies each sensed cardiac complex as either being a ventricular tachycardia complex or a supraventricular complex. The ICD 700 then determines whether the number of ventricular tachycardia complexes exceeded a predetermined threshold value. In one embodiment, when the number of ventricular tachycardia complexes exceeded a predetermined threshold value, the ICD 700 declares a ventricular tachycardia event. When a ventricular tachycardiac event is declared, the microprocessor 724 provides a signal to a cardioversion/defibrillation output circuit 744 to deliver ventricular tachycardia therapy to a heart.

In an alternative embodiment, the ICD 700 determines whether the number of supraventricular tachycardia complexes exceeded the predetermined threshold value. In one embodiment, when the number of supraventricular tachycardia complexes exceeded a predetermined threshold value, the ICD 700 declares a supraventricular tachycardia event. When a supraventricular tachycardiac event is declared, the microprocessor 724 provides a signal to a cardioversion/defibrillation output circuit 744 to deliver supraventricular tachycardia therapy to a heart.

Power to operate the ICD 700 is supplied by a battery 748. Memory 750 is also provided in the ICD 700, and is connected with the microprocessor 724. The ICD 700 further includes a transmitter/receiver 754, which can be used to communicate with the microprocessor 724 through a programmer 760 as is known.

The embodiments provided herein are intended to demonstrate only some of the embodiments of the present system. Other embodiments utilizing the present subject matter are can be appreciated by those skilled in the art. For example, the concepts of the present subject matter are expressly described in terms of cardiac complexes sensed for the QRS-wave of the heart, however, applications to other cardiac complexes, including P-wave complexes or a combination of QRS-wave and P-wave complexes, can be readily appreciated by those skilled in the art without departing from the present invention.

Also, a dual chamber implantable cardiac defibrillator can be used to distinguish SVT events from SVT events based on sensed cardiac signals. In one embodiment, the dual chamber implantable cardiac defibrillator includes an ENDOTAK single body lead catheter implanted in the ventricular region of the heart and an atrial catheter implanted in a supraventricular region of the heart. This embodiment allows for ventricular near-field signals and ventricular far-field signals, along with atrial near-field signals to be sensed and analyzed by the implantable cardiac defibrillator.

Other cardiac defibrillator systems and catheter configurations may also be used without departing from the present system. In addition to ICD systems, the present system may be utilized in external defibrillation systems and in external cardiac monitoring systems. In addition to employing endocardial leads, the present system can also utilize body surface leads.

Additionally, even though ventricular tachycardia events were discussed herein, other arrhythmic events can also be analyzed to determine the nature, or origin, of the cardiac arrhythmia using the teachings provided herein, and therefore, the express teachings of this disclosure are not intended in an exclusive or limiting sense. 

What is claimed is:
 1. A method comprising: sampling a normal sinus rhythm (NSR) QRS-complex from a NSR cardiac cycle; locating a plurality of feature points on the NSR QRS-complex based on morphological features of the NSR QRS-complex; determining a NSR template for plurality of the NSR cardiac cycles, wherein the NSR template includes a median value for each of the plurality of feature points; producing a NSR filter output from a numerical convolution of the NSR template and the plurality of feature points for each of the plurality of the NSR cardiac cycles; for the plurality of the NSR cardiac cycles, determining a median NSR filter output template, where the median NSR filter output template has a median value for each value in the NSR filter output; detecting a tachycardia event sampling a QRS-complex from a tachycardia complex during the tachycardia event; locating the plurality of feature points on the QRS-complex from the tachycardia complex based on morphological features of the QRS-complex; producing a tach-ycardia complex output from a numerical convolution of the NSR template with the plurality of feature points on the QRS-complex from the tachycardia complex; and summing a numerical difference between the values of the tachycardia complex output and the median NSR filter output template.
 2. The method of claim 1, including classifying the tachycardiac complex as a ventricular tachycardia complex when the sum of the difference between the values of the tachycardia complex output and the median NSR filter output template is greater than or equal to a predetermined sum of residual threshold value.
 3. The method of claim 2, including determining an absolute value of the difference between the values of the tachycardia complex output and the median NSR filter output template; and classifying the cardiac complex as a supraventricular tachycardia complex when the absolute value of the difference between the values of the tachycardia complex output and the median NSR filter output template is less than the predetermined sum of residual threshold value.
 4. The method of claim 3, including sampling a plurality of tachycardiac complexes, classifying each of the plurality of tachycardiac complexes as either a ventricular tachycardia complex or a supraventricular tachycardia complex, and determining whether the number of ventricular tachycardia complexes exceeded a predetermined threshold value.
 5. The method of claim 4, including declaring a ventricular tachycardia event when the number of ventricular tachycardia complexes exceeded the predetermined threshold value.
 6. The method of claim 5, including providing a signal to a cardioversion/defibrillation output circuit and delivering a ventricular tachycardia therapy from the output circuit to a heart when the ventricular tachycardia event is declared.
 7. The method of claim 4, including declaring a supraventricular tachycardia event when the number of supraventricular tachycardia complexes exceeded the predetermined threshold value.
 8. The method of claim 7, including providing a signal to a cardioversion/defibrillation output circuit and delivering a supraventricular tachycardia therapy from the output circuit to a heart when the supraventricular tachycardia event is declared.
 9. The method of claim 1, wherein summing a numerical difference includes summing an absolute value of the difference between the values of the tachycardia complex output and the median NSR filter output template. 